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Abstract 

We set up and analyze a model of radiation damping within the framework of 
continuum mechanics, inspired by a model of post-Newtonian hydrodynam- 
ics due to Blanchet, Damour and Schafer. In order to simplify the problem 
as much as possible we replace the gravitational field by the electromagnetic 
field and the fluid by kinetic theory. We prove that the resulting system has 
a well-posed Cauchy problem globally in time for general initial data and 
that in all solutions the fields decay to zero at late times. In particular, this 
means that the model is free from the runaway solutions which frequently 
occur in descriptions of radiation reaction. 

1 Introduction and main results 

The Vlasov-Poisson system is a well-known description of collisionless particles 
which interact via a field which they generate collectively. It can be applied in 
the case of particles interacting through the electromagnetic field (plasma physics 
case) or the gravitational field (stellar dynamics case). The equations modelling 
the two cases are only distinguished by a difference of sign. This description is non- 
relativistic and is only appropriate for physical situations where the velocities of the 
particles are small compared to the velocity of light. When it is replaced by a fully 
relativistic model the two cases diverge drastically. In the electromagnetic case the 
appropriate system of equations is the (relativistic) Vlasov-Maxwell system while 
in the gravitational case it is the Vlasov-Einstein system, which is much more 
complicated. 
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In classical electrodynamics it is well known that accelerated charged parti- 
cles radiate and that this leads to an effect on the motion of the particles known 
as radiation reaction. This typically leads to damping, i.e. to loss of energy by the 
particles. A similar but more complicated effect occurs in the case of the gravi- 
tational field. It is, however, hard to formulate exactly due to difficulties such as 
the nonlinearity and coordinate dependence of the equations used. There is a large 
literature concerning effective equations in electrodynamics which incorporate ra- 
diation damping without providing a full relativistic description of the field and 
sources. These effective equations usually have undesirable solutions which tend 
to infinity exponentially fast, the so-called "runaway solutions". It has recently 
been observed that nevertheless, in some of these models, the physically relevant 
solutions of the effective equation constitute a center-like manifold in phase space, 
restricted to which the dynamics is completely well-behaved. Moreover, the effec- 
tive equation is a good approximation of the full system; cf. Jl4|, [l5|| . 

In the case of the gravitational field, radiation damping is a subject of partic- 
ular interest at the moment due to the fact that gravitational wave detectors will 
soon be ready to go into operation and it is important for their effective function- 
ing that the sources of gravitational waves be understood well. (For background on 
gravitational wave detection see for instance (^] and references therein.) The most 
promising type of source at the moment is a strongly self-gravitating system of two 
stars rotating about their common centre of mass which lose energy by (gravita- 
tional) radiation damping and eventually coalesce. As has already been indicated, 
it is hard to describe this within the full theory and hence effective equations like 
those known in electrodynamics are important. Very little is understood about 
this in terms of rigorous mathematics at this time. The aim of this paper is to 
take a first step towards bringing this subject into the domain where models can 
be defined in a mathematically precise way and theorems proved about them. 

The model we will discuss has the following characteristics. It clearly exhibits 
the phenomenon of radiation damping. It does not suffer from pathologies such 
as runaway solutions. It is simple enough so that we can prove theorems about 
the global behaviour of the general solution. The particular model was chosen 
with the aim of obtaining this combination of properties. It is inspired by a model 
of Blanchet, Damour and Schafer Q for a perfect fluid with radiation damping. 
The phenomenon of radiation damping is intimately connected with the long time 
asymptotics of the system. Thus, in order to capture it mathematically, we need 
at least a global existence theorem. This seems hopeless for a fluid, due to the 
formation of shocks, and so we replace it by collisionless matter (Euler replaced 
by Vlasov) . The latter is known to have good global existence properties [jlj], [T^, 
In). Although the original motivation came from the gravitational case, the 
electromagnetic case is much simpler. Thus we use a model motivated by the 
electromagnetic case here, hoping to return to the more complicated gravitational 
case at a later date. We are not aware that the model used here has a direct 
physical application. 
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The model to be studied is defined as follows. There are two species of parti- 
cles of opposite charges, say ions ("+") and electrons ("—"). In the case of the 
Vlasov-Poisson system the motion of the individual particles is governed by the 
characteristic systems 

X+ = V + , V + =VU(t,X + ), (1) 
X-=V-, V-=-VU(t,X~), (2) 

where U — U(t,x) is the (electric) potential. The requirement that the particle 
densities / ± = f^(t,x,v) be constant along the characteristics leads to the Vlasov 
equations 

d t f + + v-V x f + +VU-V v f+ = 0, (3) 
d t r+v-V x f--VU-V v f- =0, (4) 

with iglR, a; SIR 3 , and veIR 3 denoting time, position, and velocity variable, re- 
spectively. The potential U derives from the Poisson equation 

MJ = A-K P = ATi{p + -p-), lim U(t,x) = 0, (5) 

\x\ — >oo 

where 

p ± (t,x)^J ' f ± {t,x,v)dv. (6) 

Supplied with suitable data f ± {t = 0) = (0)-(3 ) constitutes the Vlasov-Poisson 
system for two species of opposite charges; see [9, |2^] for general information on 
Vlasov-Poisson and related models. 

In order to introduce a damping effect due to radiation into (^)-@, we 
modify the characteristic equations by introducing a small additional term. Let 



D(t) = J xp(t,x)dx = J J x(f (t,x,v) — f (t,x,v))dxdv (7) 

denote the corresponding dipole moment, and replace ([!]), (0) by 

X+ = V+, V + = WU(t,X+) + eD(t), (8) 

X- = V~, V- = -VU(t,X-)-sD(t), (9) 

with an £ > small. This is to be thought of as an approximation to the full Vlasov- 
Maxwell system. It includes the electric dipole radiation which is supposed to give 
the leading contribution to the radiation reaction, cf. Jl^, p. 784]. 

The third time derivative in these equations can lead to pathological be- 
haviour and so we will modify the model by formally small corrections so as to 
eliminate it. Here we follow the procedure of H which was used to tackle the fifth 
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time derivatives which occur in the analogous gravitational problem. To reduce 
the order of derivatives on D (t), we utilize the transformations 

V + = V + -eD{t) and V~ = V~ +eb(t). (10) 

Then (|), (§) read 

X+ = V+ + eD(t), V+=VU(t,X+), (11) 

X~=V--eb{t), V~=-VU(t,X-), (12) 

where the tilde has been omitted for simplicity. The corresponding Vlasov equa- 
tions are then 

d t f + + {v + eb{t))-V x f + + VU-V v f+ = 0, (13) 

d t f- + (v-eD{t))-V x f--VU-V v f- =0. (14) 

Next we derive an approximation D^(t) to D(t). By definition of D(t) in (0) we 
formally calculate 



D{t) = J J x{d t f + ~d t f-)dxdv 

= 0(e)+ J J x(-v-\7 x f+~\7U-\7 v f + + v-\7 x f- -\7U-\7 v f-^dxdv 
= 0(e)+ ( ( v(f+-f-)dxdv. 



Thus 

D{t) = 0(e)+ J J v(d t f+-d t f-)dxdv 

= 0(e) + J J v(-v\7 x f + -VU-V v f + + v-V x r-VU-V v r^dxdv 
= 0(e)-J J vVU-(V v f + +V v f-)dxdv = 0(e) + J j "vU{f + + f~)dxdv. 
Hence we are led to define the approximation 

D [2] (t)= j J\7U(t,x)(f + {t,x.v) + f-(t,x,v))dxdv 

\7U{t,x)(p + (t,x)+p-(t,x))dx, (15) 

and we replace (|ll]), (|l^) by 

X+ = V + +eD [2] {t), V + = \7U{t,X+), (16) 
X~ = V-~eD [2] (t), V~ = -VU(t,X-), (17) 
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with corresponding Vlasov equations 

d t f + + (v + eD^(t))-W x f + + WU -W v f+ = 0, (18) 

d t f- + {v-eD^(t))-V x f- - VU -V v f- = 0, (19) 

and U is determined by (H). 

We call the system consisting of ©, ©, (§), (|), and (|l|) the Vlasov- 



Poisson system with damping (VPD), and we propose it as a model to study the 
damping effect due to radiation. We add some more comments. 

Remark 1 (a) To model radiation as we have done it, it is necessary to consider 
at least two species with different charge to mass ratios. Here we make the simplest 
choice of equal masses and two charges which are equal in magnitude and opposite 
in sign. If the charge to mass ratios were equal then the rate of change of the 
dipole moment would be proportional to the linear momentum of the system and, 
by conservation of momentum, the radiation reaction force would vanish. This is a 
well-known fact (absence of bremsstrahlung for identical particles), cf. |t], p. 411] 
or (2^, p. 201], and can also be seen from the corresponding effective equations for 
radiation reaction, cf. [|l6[ eq. after (1.9)]. 

(b) In the context of general relativity, one has to use the quadrupole moment 



(t)= I (x l x J -^\x\ 2 S i ^p(t,x)dx 



instead of the dipole moment D(t) and it is the fifth time derivative which occurs 
instead of the third before reduction 0. This leads to considerable complications. 

(c) Notice that D^(t) = for e.g. spherically symmetric solutions, whence there 
is no radiation damping in this case. 

It is the purpose of this paper to analyze rigorously long-time properties of 
classical solutions to (VPD). Therefore we first have to deal with the question 
of global existence of solutions, e.g. for smooth data functions / ± (t — 0) = of 
compact support, i.e. such that 

/■(feqnir^xIR 3 ), f±>0, and v) = for |s| > r or |v| > r 0) (20) 

with some fixed ro > 0. Since global existence is a quite non-trivial issue for Vlasov- 
Poisson like systems, we provide a complete existence proof for (VPD) in section 
[|, deriving estimates on higher velocity moments of f^ along the lines of 1 17 . This 



approach has been successfully applied to other related problems as well; cf. ^ . 
In this manner we obtain 

Theorem 1 If f^ satisfy then there is a unique solution 

/ ± GC 1 ([0,c5o[xIR 3 x]R 3 ) 

of (VPD) with data / ± (i = 0) = / ± . 
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Having ensured that suitable solutions do exist, we now turn to the decay 
estimates for quantities related to (VPD). We define the total energy 



s(t)=evn(t)+e pot (t), (2i) 



with 



£kin(t) = \J /M 2 (/ + (M,v)+/ (t,x,v))dxdv, and 

£ P ot(t) = ^J\VU(t,x)\ 2 dx, (22) 

denoting kinetic and potential energy, respectively. 
Theorem 2 Assume satisfy (HJ. Then 

£{t) = -e\DW{t)\\ (23) 
Moreover, the following estimates hold for t G [0,oo[. 

faj I|P ± WIU<C( 1 +*)"^ /orpe[l,|]/ 
0>) l|V^(t)|| J ,. a <C(l+t)- S ^ 2 /orpe[2,f]; 

( C j \DW(t)\<c(i+ty?. 

In particular, (VPD) does not admit nontrivial static solutions, and the kinetic 
energy satisfies fkin(i) £oc o,s i^oo for some Eoa >0. Moreover, if £(0)>0 and 
e>0 is small enough, then £oo > 0. 

See Section || for the proof. We note that in theorem || a slow dissipation 
of energy takes place due to the "damping term" D^(t), as can be seen from 
equation (|25|). 

Remark 2 As an aside, we include a comment on a relation to the usual Vlasov- 
Poisson system. We start with the characteristic equations (|l6|), (|l7|), i.e. 

X+ = V + +eD [2] {t), V + =VU(t,X+), 
X~ = V--eD [2] (t), V- = -VU(t,X-). 

Define 

X+ = X+, V+ = V++eD^(t), 
X~ = X~, V~ = V~ -eD [2] (t). 
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Then 

X + = V + , V + = VU{t 1 X + )+eb [2 \t)=VW{t 1 X + ) 1 
X~ = V~, V~ = -VU(t,X-) - eb [2] (t) = -VW(t,X~), 

where 

W(t,x) = U(t,x)+eD [2] (t)-x. 

Also AW — AU. Thus we obtain a solution of the Vlasov-Poisson system where 
the potential W does not satisfy the usual boundary conditions. This is similar 
to the cosmological solutions of the Vlasov-Poisson system constructed in p3| . 
They are obtained directly as solutions of a transformed system but are in the end 
solutions of the Vlasov-Poisson system with unconventional boundary conditions. 
This reformulation gives a simple way of seeing the volume preserving property of 
the flow for (VPD), since we know it for Vlasov-Poisson. It is, however, not hard 
to see it directly. 



Notation. Throughout the paper, C denotes a general constant which may change 
from line to line and which only depends on . If we consider a solution on a fixed 
time interval [0,T], and if C additionally depends on T, this is indicated by Ct- 
The usual L p -norm of a function tp = ip(t,x) over xeIR 3 is denoted by || <£>(£) Hp.^, 
and ]£<p = <p(t,x,v) and the integrals are to be extended over (x,v) € IR 3 x IR 3 , then 
we write H^WIIp-z,,- To simplify notation, an integral J always means 

Acknowledgements. We wish to thank Thibault Damour, Gerhard Rein, Ger- 
hard Schafer and Herbert Spohn for discussions and helpful advice. MK acknowl- 
edges support through a Hciscnbcrg fellowship of DFG. 



2 Proof of theorem |2j 

We split the proof into several subsections. 



2.1 Energy dissipation 

We verify (p3| ) and calculate the change of the total energy £{t) from (|2l|). Due to 
and (|19|) we have 

4in(t) = \j J V 2 (d t f + + d t f-)dxdv 

= \j J v 2 ( K -[v + eDM{t)]-V x f+~VU-V v f + 

- [v - eD [2] (i)] • Vx/" + VZ7- V„/~) dxdv 
{vVU){f + - f~)dxdv= I VU -jdx, (24) 
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where 

j(t,x)=j + (t,x)-j-(t,x), j ± (t,x) = J ' vf ± (t,x,v)dv, 

is the current. The evaluation of £ po t(t) is a little more tedious, and for this purpose 
we will use 

W*) = -^ / ^|Vt/| 2 = -^ J dx(AU)U = ^ J dx(p--p+)U, (25) 

and moreover the representations of the electric field by means of Coulomb poten- 
tials 

U(t,x) = -J -^—(p+(t,y)-p-(t,y)), 

E(t,x) :=VU(t,x) = - J dyV x j-^(p+(t,y)-p-(t,y)) 

= Jdy^t(p + (t,y)-p-(t,y)). (26) 

Then we obtain through the change of variables x <-» y and v^>w 
£pot(i) 

= \J j dxdv{(d t p--d t p+)U+(p--p+)(d t U)} 

= ~J J J J dxdydvdw {(<%/" (t,x,v) - d t ,f + (t,x,v)) 

x i^^i (/ + (*, y , w) - r (t , V , w)) 

+ (f-(t,x,v)-f+(t,x,v)) 

xj^-^(dt.f + (t,y,w)-d t r(t,y 7 w))} 

= - J J J J dxdydvdw(dtf~(t,x,v)-d t f + (t,x,v)) 

x (/+ {t,y,w)-f-(t,y,w)) 

\x-y\ 

= J J J J^^dvdw[[v-eD^(t)}-V x f-(t,x,v)-VU(t,x)-V v r(t,x,v) 

-[v + eD^(t)}-W x f + (t,x,v)-WU(t,x)-W v f+(t,x,v)} 
^{f + (t,y,w)-f-(t,y,w)) 

= J J J-^- l dv{[v-sD^(t)}-V x f-(t,x 1 v)-[v + eD^(t)}-V x f+(t,x 1 v)} 



x(p + (t,y)-p-(t,y)) 
dxdydv (Vxr^—) ■ f [v-eDW(t)]f-(t,x,v) 

V \x-y\J >- 

- [« + £^ 2] (*)]/+ (t, s, «) } (p+ (t, y) - p- (t,y)) 

dxdvWU(t,x)-{[v~eD^\t)}f-(t,x,v)-[v + eD^(t)}f+(t,x,v)} 
dxdv(VU-v)(f--f + ) + J J dxdvVU- (-eD® (t)/~ - erf® (t)f + ) 
= -J VU-jdx-e\D^(t)\ 2 , 

recall the definition of D^(t) from @. Comb ining this with we see that 

11 holds. 



2.2 Decay of the potential energy 

Here we show a t -1 -decay of the potential energy £ po t{t) from (p5|). The result 
is similar to ^] , but the proof requires appropriate modifications due to the 
presence of the term D [2] {t). 

Lemma 1 We have 

£ P ot(*) <<?(!+ i) -1 and [ f \x-vt) 2 (f + + f~)dxdv <Ct, ie[0,oo[, 



the constants being independent of £6 [0,1]. 
Proof : Denote 



R(t) = J J(x-vt) 2 (f + +f-)dxdv and g{t) = ^- J \VU\ 2 dx = 2t 2 £ pot (t). 

(27) 

Then a short calculation reveals 

R(t) = - 2t J (x-VU)pdx + 2t 2 J VU-jdx 

+ 2eD^(t)-(^J j \x-tv){f+ -f-)dxdvj. (28) 

Inserting (26) for VZ7 and with x- (x — y) |cc — j/| 3 = \x — +y- (x — y) \x — y\~ 3 , 
we see that 

y (x-VU)pdx = ~ J Updx = -^- J \VU\ 2 dx. (29) 



On the other hand, (|24|) and the energy identity (|23j) imply 
j VU-jdx = S kin (t) = -£ pot (t)-e\D^(t)\ 2 = -^j t J \VU\ 2 dx-e\D^(t)\' 



(30) 



Using (g9j) and fl3C|) in (g8|), it follows that 
t 



R(t) = - 



An 



YU\ 2 dx + 2t 2 [ - -L^- / |VE/| 2 <fe-£|£> [21 (t)f 



(x-tv)(f + -f~)dxdv 



By means of g from (|27|), this may be rewritten as 
d 



dt r K R(t)+ g (t))= g -f-2et^(t)\' 
+ 2eD^ 2 \t)- 



{x-tv)(f+ -f~)dxdv 



(31) 



Compared to (lj, p. 1412], it is now necessary to see how the two terms with 
D^(t) contribute. 

First we consider the case that t > is such that 



t 2 \D^(t)\< 



{x-tv)(f+ -f~)dxdv 



Then we obtain from (plh that 



±(R(t)+g(t)) < 9 -f +2£ \ D m (t)l 



< 9 -^ + 2st- 2 



(x-tv)(f+-f-)dxdv 

2 

(x-tv){f+ -f~)dxdv 



(32) 



However, if 



t 2 \D^(t)\> 



then (|31|) yields 



(x-tv)(f + - f~)dxdv 

g(t) 



hence ( |32| ) is verified for all i > 0. In order to have bounds below independent of, 
say, ee [0,1], we modify (|32]) to 



| (fl(()+9(())£ £ffl + ^ 



{x-tv){f + - Ddxdv 



(33) 
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To further exploit this, we next note that due to Holder's inequality and by lemma 
below with p = 



< 



(x~tv)(f + -f-)dxdv 
{x-tvf[f + + f-]dxdv 



[f+ + f-}dxdv) <CR(t), (34) 



thus by (ft 



±(R(t)+g(t))<^ + Ct- 2 R(t), *>0. 
Integrating over t£ [1,T], we see that 

R(T)<R(T)+g(T)<C+ ^j-dt + C r 2 R{t)dt, T>1. 



Therefore 



R{T)<c(l + J ^Y d ^ T ^ 1 ' 
by Gronwall's lemma. Using this in (pq), we find 



<l(T) C C+ [ T ^ldt + C C I 1 



/i * 

<c+£ g -^dt +Cj 



t gMds \d, 

i s 



f 2 



1 1 \ $(*) 



i T / t 



-dt, 



(35) 



(36) 



and consequently 

g(T)<CT, T>1, 

again by Gronwall's lemma. According to the definition of g, this proves the t^ 1 - 
decay of £ po t(i), and then (Q) shows that R{T) < CT holds as well. This completes 
the proof of lemma [l]. □ 



2.3 Some general estimates 

We digress now from the proof of theorem ^ and note some useful estimates that 
will also play a role later for the global existence of solutions, cf. theorem |l|. For 
p(z [0,oo[ define the velocity moments 

M±(t) = J Jlif^i^x^dxdv, and M p (t)= sup (m+ (s) + M~ (s)). (37) 
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Lemma 2 For t£ [0,oo[ we have 

WfHt)\\oo; X v<C and M p (t)<C, pG[0,2}. 

Proof: For fixed t g [0, oo \ let (X(s),V(s)) = (X(s;t,x,v),V(s;t,x,v)) denote the 
characteristics from (fig) associated with (Uq), i.e. 



X(s)\ = (V(s)+sDW(s)\ [X(t)\ = 
V(s)J \ VU(s,X(s)) )' \V{t) 



(38) 



Then ■§ s [f + (s,X(s),V(s))]=0 shows f+(t,x,v) = fo(X{0), V(0)), and hence the 
first bound follows. Concerning the second, 

M 2 (t) = 2 sup £ k in(s)<2 sup £(s)=2£(0) 

se[o,t] se[o,t] 

by (||). Also M ± (t) = M5 fc (0), as (s,t>)i^(Af(0;i,x,t>),V(0;i,a:,i;)) is a volume- 
preserving diffeomorphism of 1R 3 x H 3 , due to the fact that the right-hand side 
of the ODE in ( |38| ) has divergence div = div(^ j y) zero; see also lemma [l6| below 
and remark^. Observing that |w| p <l + |w| 2 forvglR 3 andpe[0,2] completes the 
proof. □ 



Lemma 3 Let f = f(x,v) S L°°(]R 3 x IR 3 ) be a nonnegative function such that 
1 1 \ v\ p f(x,v)dxdv < oo for some p£ [0,oo[, and define <fi(x) = J f(x,v)dv, x<E IR 3 . 
Then 

3 

U\\i+E, x <C\\f\\^ xv ^J J\v\Pf(x,v)dxd v y +P . (39) 
Here C depends only on p. 

Proof : The argument is well-known, but indicated for completeness. We split 
4>(x)< [ f(x,v)dv+ [ f(x,v)dv<^-^\\f\\ 00 . m +R-' [ \v\Pf(x,v)dv 

J\v\<R J\v\>R 6 

and optimize in R to find 

^)<C\\f\\^ x Jj \v\*f(x,v)dv 
whence integration w.r.t. x yields (|39"|). □ 
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Lemma 4 We have 



1 



<C\\p\\ 



p\ X ' 



3q 



In addition, 



*divT 



<c||r| 



q; x ' 



ge]l,oo[, 



Q ; x 



for smooth and compactly supported vector fields T :JR — > IR . 

Proof : The first estimate is a consequence of the classical Hardy-Littlewood- 
Sobolev inequality; see Thm. 4.5.3]. Concerning the second, we note that 
integration by parts reveals 



x-y 



x- 



47T 

divr(y) dy = — T(x 



-y|3 ----- 3 



lim 



T(x-y)-g(y)dy, (40) 

\x— y | >e 

with g(y) = j^ s G(y), where G(y) = (-ld) + j^(y(g>y)eM 3x3 . Since G is bounded 
in IR 3 \{0}, homo geneous of degree zero, and satisfies J,,_ 1 G(y)d 2 y = 0, the Cal- 
deron-Zygmund inequality III Thm. 4.31] implies that the second term on the 
right-hand side of ( f4(i| ) defines a bounded operator L q (Si 3 ) — > L q (JR 3 ); in view of 
the compact support of the T's it is not necessary that G has compact support. □ 



Lemma 5 For t€ [0,oo[ we have 

||p ± WII P;x <CM 3(p _ 1) (i)i pG[l,oo[, (41) 



as well as 



\\VU(t)\\ q . x <C\\p(t)\\^. x <CM^(t) 3 -£, ge]|,oo[. (42) 

Moreover, 

\D^(t)\<\\VU(t)\\ y Jp(t)\\^ Pe[l,oo]. (43) 
Proof : According to lemma || and lemma |^, 

ll/^(*)ll^ ;x <C'll/ ± (*)ll2S„M Q (t)'fe<CM a (t)'fe 

for all a > 0, hence ( ^l| ) holds. Due to (|26|) w e see lemma || applies to yield, for 



gG]|,oo[ and with p=^^ : together with ( |4l| ) 



l|VCr(t)|| g . s <c(||p+(t)|| 1> . x + ||p-(t)|| 1>ia! )<CM3 (p _ 1) (t)5 



Expressing^ through q, we arrive at (fl2"l). The estimate on |Dl 2 l (i)| is a consequence 
of ( |l5| ) and Holder's inequality. □ 



13 



2.4 Proof of theorem Q (completed) 

From lemma [j] we additionally obtain, analogously to fl2]| , the following informa- 
tion. 

Corollary 1 Under the assumptions of theorem^, we moreover have 

||p ± Wll f;x <C(l + t)- 3/5 , te[0,oo[, (44) 

and 

||V{/(i)||i5 ;2; <C(l + tr 3/5 , tG[0,oo[. (45) 
Proof : Using lemma ^ and lemma [j] we can split 

P ± {t,x) < / f ± (t,x,v)dv + R~ 2 / (x-tv) 2 f ± (t,x,v)dv 

J {v.\x-tv\<R} J {v:\x-tv\>R} 

< CR 3 t- 3 + R- 2 J {x-tv) 2 {f + + D{t,x,v)dv<CR 3 r 3 + CR- 2 t, 

and then choose the optimal R = t 4 ^ 5 to obtain (H). Concerning (ff5|), this follows 
from Q) and the first inequality in (|i|) with q = ^f. □ 

Lemma 6 Assertions (a)-(c) of theorem^ are satisfied. 

Proof: From lemma |2| and corollary [l] we know that ||p ± (t)|| 1 . :E < C, and moreover 
\\p ± (t)\\s :x <C(l+t)- 3 / 5 . We may also estimate || VU(t)\\ 2 . x = V&^£ po t(t) 1/2 < 
Cil + t)- 1 / 2 as well as || VU(t)\\ is. x < (7(1 + i)~ 3 / 5 by lemma and corollary 0. 
Hence the general interpolation estimate 

p qi 92 P qi Q2 

yields (a) and (b). For (c), we use J43| ) with p = 5/2 andp' = 5/3, (a), and (b) to 
see that 

|^ [2] (*)l < <C(l + i)- 19/35 (l + *)- 3 / 5 = C7(l + t)- 8/7 , (46) 

as was to be shown. □ 



Remark 3 The estimates derived thus far suggest that the optimal decay rate for 
D^(t) be \D^{t)\^t- 3 ' 2 rather than \D^ (t)\ ~^ 8 / 7 , for the following reason: 
by Holder's inequality and lemma [j] the quantity I(t) — J f(x — vt)(f + — f~)dxdv 
satisfies 

1 /2 

\I{t)\<c(J J(x-vt) 2 (f+ + f-)dxdv^j <C(l + t)^ 2 , 
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thus we might expect J(t)~f 1 / 2 . On the other hand, direct calculation shows 

Ht)=(e J J(f + + f-)dxdv~?)DM(t)~(-t)DW(t), 

whence we should have \D^(t)\~t~ 3 / 2 . This decay would also be obtained if it 
were possible to use theorem ||(a) and (b) with p=^f- and p = yy , respectively, 
since then (fl3j) would yield 

\D^(t)\<\\VU(t)\\^Jp(t)\\^ x 

5(15/ll)-3 3((15/4)-l) .,,„ 

<C(l+t) WD] (1+t) 2(15/4) =C(l+t) d/2 . 

However, the necessary decay estimates for such p-norms of VU{t) and p(t) could 
not be proved. 



Corollary 2 There are no nontrivial static solutions of (VPD), and £kin(i)^ 
£ao > as t—f oo. If £(0) > and e > is sufficiently small, then £oo > 0. 

Proof: If (VPD) had a static solution f ± (t) = then £ pot (t) = 0, whence VU 
0. This in turn yields D^(t) = by definition. Consequently, the Vlasov equations 
©, (|l9|) reduce to a t /± +u- V x / ± = with unique solution ^(x^) =f±(x- 
vt,v). But then we see p^{x) = J f^(x — vt,v)dv — t~ 3 f f^ z (w,t~ 1 [x — w])dw, show- 
ing as t— >oo that the solution has to be trivial. To prove the assertion concerning 
£kin(i), note that, since £(t) is decaying by (p3|), £(t)^£ oc >0 as t— >oo. But 
=£kin(i) +£ P ot(£) and £ pot (t)^0, hence the first claim follows. For the sec- 
ond, denote C\ the constant on the right-hand side of (p!6|). Since all bounds are 
derived from lemma Q, we note that C\ is independent of £ € [0, 1] . Integrating ( |23| ) 
yields 

£km(*)+£pot(*) = £(0)- £ / l^ [2l (s)| 2 ds, 

Jo 

thus as i^oo, provided that foo =0, by ( fi6| ) 

/•oo />oc 

£{0)=e \D^(s)\ 2 ds<Cfe (1 + S y 16/1 ' ds = {7Cf/9)e. 
Jo Jo 

So if we choose e < (9/7C 2 )£(0), then necessarily > 0. □ 



Taking into account Section 2T , lemma § and corollary |2| we note that the 
proof of theorem || is complete. 

Remark 4 With regard to corollary 0, £oo = limt^oofkin^) > was to be ex- 
pected, since otherwise the particle velocities would have to tend to zero. It is, 
however, not surprising that at late times, when we are in a small data regime and 
the radiation reaction force is getting small, the solution behaves like a solution of 
the Vlasov-Poisson system with small data. In that case the particles travel with 
constant non-zero velocity at late times, as shown in p|. 
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We note a further consequence of the foregoing estimates. 
Corollary 3 If£ oo >0, then 

C 1 t-C 2 <J J{x-v)(f+{t 7 x,v) + f-(t,x,v))dxdv<C 3 {l+t), fe[0,oo[, 

for constants Ci,C2,Cs>0. 

Proof : Denote S(t) = J J (x-v)(f + + f~)dxdv. In view of lemma [I] we obtain 

2t 2 £ kin (t) = J J{x~vt) 2 (f + + f-)dxdv- J J x 2 (f+ + f-)dxdv + 2tS(t) 
< C{l + t) + 2tS{t), 

whence t 2 £oc < C(l + t) + 2tS(t) for t large enough. To prove the upper bound, 
note that Q(t) = / f x 2 (f+ + f~)dxdv satisfies Q(t) = 2S(t) + 2eD^ 2 \t) ■ J J x{f+ - 
f~)dxdv, as follows by a straightforward calculation. Therefore utilizing Holder's 
inequality we obtain 

Q(t) < 2Q{t) 1 ' 2 (2£ k Ut)) 1/2 +C{l + t)- B ' 7 Q{t) 1 ' 2 < CQit) 1 ' 2 . 

Consequently, Q(t) < C(l + t) 2 , and this in turn yields, once more by Holder's 
inequality, \S(t)\ < Q(t) 1 / 2 {2£ kin {t)) 1 / 2 < C(l +t). □ 

3 Existence of solutions 

As mentioned in the introduction, the proof follows jf7j . The idea is to decompose 
the field E — VU = E\+ F in a "far field" F that is small, and in some com- 
plementary part Ei which is of higher regularity than E itself. (More precisely, 
\\Ei(t)\\ p . x <C for every p£ can be achieved.) According to this splitting, 

we write the Vlasov equations (|l8) and in the form 

d t f + + (v + eD^{t))-V x f++F-V v f+ - -E 1 -W v f+, (47) 
d t f- + {v-eDW{t))-V x f--F-V v f- = Ex-V v f-. (48) 

Since F is small, the characteristics of e.g. should behave as 

X{s)^x+{s-t)v + e [ D [2] (r)dT, V(s)w«, (49) 



which is close to a free streaming, at least in case were not present. Writing 
p^(t,x) as a suitable integral over characteristics, it then turns out that in order 
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to derive the necessary estimates for global existence (on higher moments), it is 
possible to use a rigorous form of (fly). In particular, one may verify that 



and 



dx 



■t\, 



as is important to transform away the characteristics. The main point to note here 
is that the term with drops if we take derivatives in (^) w.r.t. x or v, and 
hence the arguments from [jl7| can be expected to carry over. Having derived the 
higher moment bounds 

M m (t)<C, te [0,1*1, mG]3,^[, 

then a standard argument yields the global existence of classical solutions for 
(VPD). 

It should finally be remarked that we did not succeed in generalizing the 
proofs of global existence for the usual Vlasov-Poisson system that bound the 
increase in velocity along a characteristic; see Jh], |24|, [22fl , and [^o) for a recent 
application. The reason for this is that, when estimating the "ugly" term, an X(s) 
will appear, which in our case will lead to the expression D^(s) that could not be 
bounded well enough to make the proof work. 



3.1 Local existence 

Similar to the case of the usual Vlasov-Poisson system, cf. jjj , where this is con- 
tained implicitly, an iteration scheme may be set up to yield the local existence of 
a solution and a criterion when a local solution in fact will be global. 



Theorem 3 Suppose /q satisfy (2C). Then there exist unique solutions 



/ ± GC 1 ([0,T,[xIR 3 xIR 3 ) 

of M), and ( fiH j with data / ± (i=0) = , on a maximal time interval 

of existence [0,T*[. If moreover 



P = P H 



P~ 



with P ± (t)=sup||i;| : 3a; e H 3 3s € [0,i] : (x,v) € supp/ ± (s)|, 

(50) 



is bounded on [0,T*[, then T^~oo. 

We will not go into the proof of this result. 



3.2 Some preliminary estimates 

We first need to derive some a priori bounds. For this we consider a classical 
solution of the system that exists for times t£ [0,T]. Note that all estimates from 
the previous sections remain valid on any interval where the solution exists. 
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Lemma 7 For t£ [0,T] 



M p (t)<C T (l+ sup ||Vtf(s)fc V pG[l,oo[. 
v «e[o,ti - 



In addition, 

MpUXCr^l + M^+^^t)^^, pe [l,oo[. 
Proof: Recalling j37|), from ( fL8| ) and Holder's inequality it follows that 
j t M+(t)= J J\v\p[~[v + eD^(t)}-V x f+-VU-V v f+}dxdv 



(51) 



JP-2, 



v-VU)f + dxdv 



< C\\VU(t)\\ 3+p . a 



\v\v- l f + {t,;v)dv 



3 + p . 

2+p ' 



Now 



<CM±(ty+* 



by an argument similar to the proof of lemma [|, and this yields 
d 



dt 



< 



C\\VU(t)\\ 3+p . x M p {t)*& . 



Since M p {-) is increasing, it is differentiable a.e. in t, with 



dt 



M p (t)< sup 
se[o,t] 



dt 



M+(s)+M p (s) 



2+P 



<C sup \\VU(s)\\ 3+p . x M p (t)^. (52) 

s£[0,t] 



Integration of this differential inequality gives the claim. Finally, for (^l|) we ob- 
serve that by the first part and ( (42| ) with q = 2>+p 

M p (t)<C T + C T sup \\VU(s)\\lXl. x <C T + C T Me^o{t)^r {3+p \ 

«G[Q,t] 3+5 



and 6£_9 =3 (3+2p). 

3+9 v 6+p / 



□ 



3.3 Estimates for higher moments 

For R>0 choose a radially symmetric function xr & Cq°(JR 3 ) with Xfl(^)G[0,l] 
for irelR 3 , xr{ x ) = 1 Ior \ X \<R, an d = f° r |a:|>2.ft. Correspondingly we 
decompose the electric field E(t,x) from ( |2q ) as 



£(i,a:) = .E 1 (i,aO + .F(<; ) a:), 



18 



with 



Ei(t,x) = I x (x-y)f-^ip+(t,y)- p -{t t y))dy 
\x-y\ a 

= -[xV±)*(p + (t)-p-(t))(x). (53) 



Some useful estimates on E\ and F are stated in lemma [T^ below. Then we 
write the Vlasov equations (|l8|) and (|l9|) in the form (|^) and (|||). This can 
be used to derive a representation formula for p^(t,x), and for simplicity we will 
consider only p + (t,x). We fix x,vefR 3 and te[0,T], and denote (X(s),V(s)) = 
(X(s;x,v),V(s;x,v)) for s€ [0,t] the solution of the characteristic system 



X(s)\ = (-V(s)-eDW(t-8)\ [X(0)\ 
V(s)J \ -F(t-s,X(s)) )' {V(0)J 



(54) 



associated with (^) . Since 
^-[f+(t-s,X(s),V(s))} 

OS 

= -d t f + (t- s ,X(s) 1 V(s))^[V(s)+eD^(t-s)}-d x f + (t-.s,X(s),V(s)) 

~F(t~s,X( S ))-d v f + (t-s,X(s),V(s)) 
= E 1 (t-s,X(s))-d v f + (t-s,X(s),V(s)) 

by (|47]), it follows through integrating f*ds(...) and f dv(...) that 
p+(t,x) = j dvf+(X(t),V(t)) 

r ds J dvE 1 (t~s 1 X(s))-d v f + (t-s 1 X(s),V(s)) 
dvf+(X(t),V(t))~ [ ds j dvdbr v [E l f + ]{t-s,X{s),V{s)), (55) 



where /+ = /+(t=0), and [E 1 f+]{T,X,V) = E 1 (T,X)f+{T,X,V); note that the de- 
pendence on x and v in ( p5| ) enters via X(s) and V"(s). To rewrite ( |55] ) appropri- 
ately, define 

G(X,F) = [.E 1 / + ](i-s,J£,\/) and G(^,u) = G(X(s;x,i;),V r (s;a;,u)). 

By lemma |l6| below we then have G(X,V) = G(x(s;X,V),v(s;X,V)), and conse- 
quently 

3 
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where div x (|^)- (5 = ^^=1^ gtr(§^)- Utilizing this in @ and integrating by 
parts w.r.t. v yields 



P + (t,x) = J dvf+(X(t),V(t))-div x J dsjdv^ 



dx 
dV 



■G 



ds / dv 



dx 



divJ — -G + divJ — -G 



=: </>+(t,x)-div K r+(t,a;) + i? + (t,iz;). 



dv 



dV 



Similarly, we have 



p (t,x) = (/> (t,x)-div x T (t,x)+R (t,x), 



(56) 



(57) 



with the corresponding functions <p , T , and R 

Next we derive some estimates on (j>^ , T + , and R + . 

Lemma 8 For t 6 [0,T] we ftave 

II^WIIi+s ;:c <Ct (m>0) and H^WUa^).* < C T (m>3). 

Proof : We can apply corollary |?] below with s — t and r = to obtain the first 
bound. Concerning the second, note that m>3 implies 3(|i^)<^i™. Whence 
it suffices to bound the support of £Ei— > </>q (t,a;) = J dv (X(t;x,v),V(t;x,v)). To 
do so, recall from ( |20| ) that f£ (x,v)=Q for |a:|>ro or |u|>ro- From the proof 
of corollary |?j we know \V(t) — v\<Ci, C\ depending only on T. Thus J^\X(s) — 
x\ < \V(s)+eD^(t-s)\<C(l + \v\) by (§4|) and theorem |(c), whence \X(t)-x\ < 
C 2 (l + \v\). Then, if \x\ > G 2 (l + C ± + r ) +r =: n and \V(t)\<r , we have |w| < 
\V(t)~v\ + \V(t)\<d+r and therefore \X(t)\>\x\-\X(t)-x\>r . This yields 
f£(X(t),V(i))=0, and thus 0^(i,as)=O for \x\>n andte[0,T]. □ 

Next wc turn to bound T + (t,x) = J^ds J dv(Sy -G). 
Lemma 9 For t € [0,T] and any to G]0,T] we /iawe 

||r+(t)|| 3+m;x <G T t ^(l + M m (t) (6—)%+™)) +Or(l + |Into|)(l + M ro (t)»fe 
wit/i 3 < m < yj . ifere Gr does not depend on to . Moreover, 

\\r + (t)\\ 3+m . x < c T t^ (i+M m (t) (6-*.)%+-) ) , iG [ ,T]. 

Proof : We first note that 



(58) 



\T+(t,x)\<C T [ dsi 
Jo 



dv 



[Bi/+](t-s,X(«;a; ) i;),V(«;a;,t;)) 



(59) 
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in view of the second estimate in (^) below. Next we observe that due to lemma 
and the first estimate in (p7[) we have, with - + ^j = l, 



dv 



[3i/+](t-s,X(s),V(s)) 

d«|Bi(t-s,X(s))| 



< sup ||/+(r 

tG[0,T] 



< C T s~ 



dXlExit-s.X)^ 



<C T s~- sup WE^W 

re[0,T] 



dvf+(t-s,X(s),V(s)) 
dvf + (t-s,X(s),V(s)) 
dvf+(t-s,X(s),V(s)) 



(60) 



In fl59| ) we then split f*ds = f Q °ds + J^ Q ds with i o G]0,i]. 

To bound the first term, we choose r' = 6 — m < 3, whence r = §55 > § ■ Then 
with p= 6 "~ 9 > 1 we find ^±w = ^t £ . Moreover, Kr<^ by the choice of m. 
Hence by lemma |15| and corollary |7] 



<iss / dv 



3+m; a; 



T £ 



sup ||£i(r)|| 

re[0,T] 



x sup 

r£[0.f] 



< C7 T t r l + Af p (t)»+» 



ri;i; 



d«/+(t-T,JC(r),y(r)) 



note that Jensen's inequality has been used for the first estimate. Utilizing the 
relation 3(^±^) = m and (51), we can bound M p (t) by means of M m (t), as 

M p (t) < C T (l + M m (t) = C T (l + M ro (i) ' 
Thus we have shown that 
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|8-n 



<iss / ch; 



[E 1 /+](t- S ,X( S ;-, U ),F( S ;-,v)) 



<C T to J (l + M m (i) <•-»)%+«)). 



3+m; a: 



(61) 



As far as the second part of the integral is concerned, we now make use of 
[E 1 f + ]{t-s,X{s-- 1 v),V(s--,v)) 



0|) with r' = 3 and r = |. Then 



dss I dv 



3+m; x 
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< C T (1 + | Into|)( sup WExWU.a 

V rG[0,T] 2 ' 



(62) 



x sup 

r£[0,i] 



cb / d V / + (<-r,X(r),y(r)) 



< C T (l + |lnf |)(l + M m (t)^), 



(63) 



again by lemma [l5] and corollary ^. Summarizing ( |6l| ) and (|62] ) , we see that the 
first asserted estimate holds. To verify ( ^8| ) it sufficient to follow the argument 
just elaborated and to note that for to = t the contribution of the J t ds(. . .)-part 



of ( |59| ) drops out, whence we simply use ( |61| ) for to = t. 
Finally we need to consider 



□ 



R+(t,x) = J ds J dv 

in @. 

Lemma 10 FortG [0,T] we have 



dx 



dV 



div x — •G+divJ — G 



dv 



dV 



\R + {t)U {m) , x <c T , 



with me[0,f]. 



Proof : Using ( p8| ) below, and (pQ ) with r'= and r = ^r> we estimate 



|i? + (t,a;)| < C* T / ds / cfo [E 1 /+](t-s,X(s;ar,?;),F(s;x ! t;)) 



<C T ( Bup r ||Bi(r)||^ ilB ) / dss--^ I dvf + (t-s,X(s),V{s)) 



<C T / dss-^[ dvf + (t-s,X(s),V(s)) 



the latter according to lemma |l5|. Hence due to corollary ^, with p determined 

through 3+p — 27 ( 3+m ) 
mrougn 3 - 13 { 6+m ), 



R + (t)\\ 3( ^ )lx < C T sup 



re[0,t] 



dx[ / dvf+(t-T,X(T),V(T)) 
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3 + my 
6 + m J 



< Or l + M^TO =C T l + M 3M (f)ww , 
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where Q = §(|±^). Since 0<3(a-l) = ^(£t^E)<2 by choice of m, the claim 
follows from lemma ||. □ 

The foregoing estimates, and analogous ones for r - , and Rr , can be put 
together to yield the following result. 

Lemma 11 For tE[0,T] andm£]3,jj[ we have 

\\VU(t)\\ 3+m . x < C T tf^ (l + M m (t) 

+C T {l + \hxt \)(l+M m {t)^), 
with fo€]0,t] being arbitrary. Here Ct does not depend on to. Moreover, 

\\VU(t)\\ 3+m . x <C T t^ (l+M m (t) (•-»)%+»)), tg[ 0l T]. (64) 
Proof : Due to ( p6| ) , (^6|) , and (|57|) we may write 

vt/(^) - - (v^ ) * (W(t) -C(<)] + (*) - R~(t)] 

+div x [T-{t)-T+(t)]y x ). 
Therefore lemma [| implies for m€ [0,oo[ that 

\\VU(t)\\ 3+m . x < c(||^(t)|| 8( ^ )iie + |^o (t)|| 3(SS);a 

+ \\R + (t)\\ 3(mhx + \\R-(t)\\ 3mhx ) 

+c(\\r + (t)\\ 3+m . x +\\T-(t)\\ 3+m]X ). 

Due to lemmas |[ ^, a nd [To| we thus obtain the first desired estimate. Concerning 
(|64|), we rather apply (|5§|) than the first estimate from lemma |. □ 

This in particular can be used to derive a short time bound on M m {t). 

Corollary 4 For me]3,yj[ there exist t\ g]0,T] and Ci>0 ('feoi/i depending on 
T ) such that 

M m (t)<d, *e[o,ti]. 



Proof: Combining (52) with (J6J) and observing is-™ <Ct yields 



^M m (t) < C sup ||VE/( S )|| 3+m . x M m (*)^ 



< C r (l + M m (t) <6- m )%+m) ) M ro (t) i 
23 



Integration of this differential inequality gives a local bound on M m (t), that, how- 
ever, fails to extend to all of [0,T] due to (6 _ m) 9 (3+m) + = g=* > 1. □ 



Corollary 5 For£€[0,T] and 77i€]3,yj[ we have 

||VC/(t)|| 3+m . x <C T (l + |l n M m (t)|)(l + M m (t)^). (65) 

Proof : Note that we may assume M m (t) > 1 for £G [0,T], since otherwise M m (£) 
simply can be replaced by M m (t) + 1. Set 

in lemma [llL with t\ from corollary |J. If £e[£i,T], then £o<£, and therefore 

(^) + 7« s = aj~ shows that (BE) holds for £e[£i,Tl. On the other 

hand, if £<E [0,£iL then M m (£)<Ci and (||) imply || VU(t)\\ 3+m . x < C 2 for some 
C2 > 0. Hence (^|) holds as well in this case if we choose Ct >C 2 - n 

Theorem 4 For te [0,T] and m€]3,fi[ we have 

M m {t)<C T . 

Proof : By (52) and due to corollary |5| we see that 



^■M m (£) <C sup ||VC/(s)|| 3+m . x M ro (£)^<C' r (l + |lnM m (£)|)fl + M m (£) X ) 
"i se[o,i] ' v ' 

< C T (l + |lnM m (£)|M m (£)). 
Integration of this differential inequality yields the claimed estimate. □ 



Corollary 6 For £€ [0,T] we have 

no 

II^(*)II P;B <C T , pe}2,-[ 



Proof: This is a consequence of ( |4l| ) and theorem^, since to = 3(p — 1) e]3, |j[ 
corresponds top e] 2, f| [. □ 
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3.4 Global existence of solutions 

We start with some preliminary (well-known) observations. Recall the definition 
of P ± (t) from (50), and also that [0,T*[ is the maximal interval of existence, 
cf. theorem ||. 

Lemma 12 We have 

P ± (t)<P ± (0)+ [ ||Vl/( a )|| ds, fe[0,T„[. 
Jo 

Proof: Assume e.g. (x,v) € supp/ + (s) for some xeIR 3 and sG[0,t]. From the 
proof of lemma H we have f + (s,x,v) = (X(0;s,x,v), V(0;s,x,v)), with (X ,V) the 
characteristics from (|38|). This means that (x,v) = (X(s;0,Xq,vq), V(s;0,xq,vq)) for 
(x ,vq) € supp/o 1 ". Hence 



V(T;0,x o ,v o )dT 



<P+(0)+ / \\VU(r)\\dr 



by the characteristic equation for V. □ 
Next we need to derive a bound on || V?7(t)|| 0o . a .. 



Lemma 13 For a > | we have 



llv^(t)IU x <c £||^)llooJ Ell^WlUx) . ^^ T *l 

where — + = 1. The constant C depends only on a 
Proof : With R > we estimate from (53) 

|VC7(*,ar)| < / ^L^L+ {t , y ) + p -(t,y))) 

J\y-x\<R \ x — y\ v ' 

^<S^ (t)lu -) fl+ (^^) s (S llp,(t)lu 



#-2 



Choosing the optimal i? yields the claim. □ 
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Lemma 14 We have 

\\p ± (t)\\^ x <CP ± (t) 3 , fe[o,r„[, 

with G depending only on the data. 

Proof : By definition of P + (t) and bounding ||/ ± (i)|| 00 . xv as in lemma 0, it follows 
that 

P + {t,x)= I f+^X^dvKCW^iM^^P+ttfKCP+itf, 
J\v\<P+(t) 

as was to be shown. □ 

Using the criterion from theorem |^ and by means of corollary ^ we are finally 
going to complete the proof of theorem |l|. 

Proof of theorem |l]: Assume T* < oo in theorem |^. All the estimates on the 
moments remain valid if [0,T] is replaced by [0,T*[, since in the constants only 
terms of the form CT r , T™, and e CT * do enter. In particular, 



|p ± (t)||ai ;x <C, tG[0,T.[, 

lary |6J, where here a 
>1 , which corre 

[jjthat 



by corollary 0, where here and below the various constants C depend on T*. Choos- 
ing a = y| > ^ , which corresponds to a' = || , we deduce from lemmas [jj], [D| and 



p + (t) < p + (o)+c / Ve ii^iu) 11 (e 11^%;: 

<P+(0) + cf ds(p+(s) 3 +p-(s) 3 ) TT . 



This implies 

P(t)<P(0)+C P(s)rtds, t€[0,T»[, 
Jo 

and hence the boundedness of P on [0,T»[. □ 

We remark that we did not try to reduce the exponent yj so as to obtain the 
optimal power of P. 

3.5 Some technical lemmas 



Lemma 15 Define Ei and F as m (\5$) . Then we have \\Ei(t)\\ p . x <C for t € [0,T] 
and pG [1, as weH as 

II^WIIoo^ + llV^WIL^ + IP 2 ^)!!^^^- 2 , te[0,T]. (66) 
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Proof: From (||) and Young's inequality @ p. 29] with 1 + 1 = 1 + 1 we obtain 



ll^l(*)llp;x = 







< 








p; a; 



\p + (t)- P -(t)\L, 



q; x 



We have x(') ]~^3 € L 9 (IR 3 ) for ge[l,f[ and ||pf±=( t )|| r;x < C for r€[l,f] due to 
theorem |^(a). Combining those values for q and r, we see that we need to have 
p€[l,^r[' The bounds in (|6^) are obtained by observing that 



F(t, X ): 



l-xix-y) 



U !l) (p+(t,y)-p-(t,y))dy 



\x-v\>r n ' \ x ~ v\ 

where x&C^°(JR 3 ), and moreover J p ± (t,y)dy = J J f ± (t : y,v)dydv <C 



□ 



Lemma 16 For fixed t €]0,T] and s G[0,t] consider the map 

Z{s) : M 6 3 (x,v) 1 ^ (X (8;x,«), V(»;a;,t;)) = (X(s),F(s)) € IR 6 , 
where (X(s),V(s)) is the solution of the characteristic system |5^[j ; i.e. 

x( S )\/-ns)-^ [2] (*-s)A ^(o) 
t>( S ); ^ -F(t- s ,x( S )) ^(0) 

TTien Z(s) is a volume-preserving diffeomorphism, and 



dct 



<c s - 3 , 



<Cs, 



as well as 



dvi\dV k { ' 



<C, l<i,j,k<3, 



(67) 



(68) 



if R>0 is chosen sufficiently large (depending on T). Here 

Z(s)- 1 : (X,V) ~ (x,v) = (x(s;X,V),v(s;X,V)) 

is the inverse of Z(s), and the constants C do depend only on T, but not on 
(t,s,x,v). 

Proof : As the right-hand side of the characteristic system has divergence div = 
div(x.v) =0, the first claim follows; cf. also remark || Moreover, we have 



<cirv, 



<CirV, 



dV 



dV 



<CR~ 2 s 2 



<CR~ 2 s, 



(69) 
(70) 
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where Id denotes the unit matrix in IR 3 . E.g. to validate the estimate on the v- 
derivatives one can introduce, following ||, the function 0(s) = ^jf-(s) + sld and 



calculate that 0(0) = 0, 0(s) = -f|(s)+ld, 0(0) = 0, as well as 

J>(s) = ^(t-s,X(s))[4>(s)-sld}. 

Here the important observation is that 
d (dX 



0( S ) 



ds V dv 



is) 



+ Id = -?-X(s) +Id = ^- f - 7(s) - eDPl (t - s)^ + Id 



dV 

at; 



since the term with e£)[ 2 l(i — s) simply drops through the ^-derivative, and thus 
the same general argument can be used as in the case without D^. Then we 
may write 0(s) = J Q S (s — T)0(r)<iT and utilize ( |66| ) to derive ( |69| ) by means of a 
Gronwall argument, whereas ( [fo| ) is obtained in the same way using the function 

4>(a) " 



- ^(s)-Id instead. 
Then 



det 



dX 
dv 



00 



det 



dX 
dv 



(s) + sld 



-Id 



together with (69) and the continuity of the map Ay—> \det(A)\ at A = Id yields 
that for s £ [0,T] and R> large enough, 



, . dX . . 
det(^( S ) 



>-s J 



and this proves the first estimate in (|67|). As what concerns the bound on |^(s)|, 
employing the chain rule it follows that 



dx 



dV 

dx 
dX 



dx 



ax 



ax 



<9v; 



9i> 



(*) = 



Id- 



9x 
97 



00 



oX 

9a; 



00 



Choosing R > sufficiently large, we find from (p9h and (J70f) that 



9t; 



ax 

9x 



whence by (|69|) and ([70| 



So 



So 



<C| l + i?" 2 s 
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and this gives |^(s)| < Cs, for R> large enough. 

Finally, the estimates on the second derivatives in (^8|) are more tedious, but 
verified in a similar way. □ 



Corollary 7 If (X(s),V(s)) is a characteristic curve, cf. lemma^jj^, and s,tG 
[0,T], then for pe]0,oo[ 



/ ± (T,X(«;-,t;),V(«;-,iO)*> 



<a 



' 1 1 / ± to 1 1 5l„ ( M o (t ) * + M p (r ) sfe ) , 



with Mp(-) from ([5 

Proof: Utilizing lemma||with f(x,v) = f ± (T,X(s;x,v),V(s;x,v)), it follows that 
the left-hand side is bounded by 



3 
3 + p 



From the characteristic equation for V(s) we obtain that -^\V(s) — v\<\\F\\ lx> . xt < 
C, whence \V(s) - v\ < CT = C T , thus \v\ p <C(l + \V(s)\ p ). Using this estimate 
and the fact that Z(s) is a volume-preserving diffeomorphism yields the claim. □ 
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